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Dynamic spin correlation functions {Sf{t)Sj) for the ID S = 1/2 XX model H = —J'Zi{SfSf+i + 
^i+i} calculated exactly for finite open chains with up to A'^ = 10000 spins. Over a certain 
time range the results are free of finite-size effects and thus represent correlation functions of an 
Ch ' infinite chain (bulk regime) or a semi-infinite chain (boundary regime). In the bulk regime, the 

long-time asymptotic decay as inferred by extrapolation is Gaussian at T = cxd, exponential at 
< T < oo, and power- law (~ t-^/^) at T = 0, in agreement with exact results. In the boundary 
' regime, a power-law decay obtains at all temperatures; the characteristic exponent is universal at 

T = (~ t~^) and at < T < oo (~ t'^^^), but is site-dependent at T = oo. In the high- 
temperature regime {T/J S> 1) and in the low-temperature regime {T/J <C 1), crossovers between 
different decay laws can be observed in {Sf{t)Sj}. Additional crossovers are found between bulk- 
. type and boundary-type decay for i = j near the boundary, and between space-like and time-like 

J — ' behavior for z / j. 

o ■ 



I. INTRODUCTION 



The long-time behavior of correlations for quantum many-body systems in general and for quantum spin systems in 
particular has been a notoriously difficult subject of theoretical research. Approximation schemes tend to have little 
' reliability in this field. There exist very few exact results for nontrivial cases, and many of them exhibit non-generic 
Q ' features for one reason or another. In classical many-body systems, the long-time correlations can be investigated by 
O means of computer simulations, but no practical quantum counterpart of that approach exists. 

^ In some cases, useful conclusions on the long-time behavior can be drawn from a moment expansion via rigorai^ 

bounds, but the time intervals over which stringent bounds can be established are often too short for that purpose.trEl 
The continued-fraction analysis based on the same number of moments can be used to predict the exponent of a 
power-law infrared singularity in the frequency domainj-|hut this approach tends to be insensitive to subtle changes 
in the long-time decay if it does not involve power laws.aiJ 
The one-dimensional (ID) 5 = 1/2 XY model, 



N-l 



HxY^-J2 {J-S^St,^ + JySfSf^, + hS^} , (1.1) 

1=1 

is one of the very few many-body systems with nontrungLLdynamics for which time-dependent correlation functions 
have been calculated exactly at-zero&lla and nonzeroE3T3 temperatures. This model is equivalent to a system of 



noninteracting lattice fermions.tj'EZI The spin correlation function {S^{t)Sj) is a simple fermion density correlation 
function, and the function {Sf{t)Sj) can be reduced to a determinant whose size increases linearly with the number 
of sites between i or j and the nearest boundary of the chain. 



The focus of this study is on the dynamics of the spin ^-components for the special case = Jy, h — of (1.1) 
- the XX model. For quite-SQiiBe time it has been known that the rfwiction {Sf{t)Sf) of the infinite system exhibits 
a Gaussian decay at T = oot2rtil and a power-law decay at T = O.ffl A more recent study states that the long-time 
asymptotic decay of the same correlation function is exponential at finite nonzero temper atures.cJ Numerical evidence 
for exponential decay was also found for {Sf{t)Sf) in the XXZ model at T = ooB 
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For a semi-infinite XX chain the available evidence indicajhes that the function {Sf{t)Sf) exhibits a power-law. 
decay at all temperatures. Rigorous results exist for T = oo,u'E3 and the result of a finite-chain study for T = 0,E£l 
but no results for < T < oo appear to have existed prior to this study. 

The purpose of this paper is (i) to fill in the missing links on the question of long-time asymptotic behavior and (ii) 
to elucidate various kinds of crossovers between the different decay laws that can be found in the autocorrelation and 
paircorrelation functions {Sf{t)Sj). The determinantal expressions for these functions have been known for a long 
time.B However, only with today's advanced computer technology can they be evaluated for systems large enough to 
yield data from which conclusions can be drawn with some confidence about the long-time asymptotics for infinite 
and semi-infinite chains at arbitrary temperatures .c2l In Sec. II we describe the method used for our analysis, in Sec. 
Ill we present our results for the infinite system, and in Sec. IV we discuss boundary effects. 



II. FERMION REPRESENTATION 

The Jordan- Wigner transformationEllli 



S^ = aW^l, (2.1) 



S+ ^ -l-^a] , S- = a,(-l)Sr.i , (2.2) 

between the component and ladder operators Sf , — Sf ± iSj for an array of localized spins with 5 = 1/2 and 
the creation and annihilation operators aj, of an array of fermions converts the Hamiltonian of an open-ended XX 
chain, 

Af-l 

Hxx = -JY. {^^^^+^ + ' (2.3) 

i=l 

into a Hamiltonian of nonintcracting fermions 

^ JV-l 

Hxx = -^Y. + al+io.) • (2.4) 
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The energies of the one-particle eigenstates are 



Ek^-Jcosk ; fc = y ; iy = l,...,N. (2.5) 



The spin correlation functions {S^{t)Sp are then in essence density correlation functions for the bandj2.5) of 
free fermions. Their characteristic long-time asymptotic behavior at zero and nonzero temperaturesaQlij is a 
consequence of the band-edge singularities in the one-particle density of state and (for T = 0) the singularity at w = 
generated by the Fermi function. For the spins at the boundary of a semi-infinite chain, different power-law decays 
of lS^{t)Sj) pertain to T = (~ if both i and j are odd and ~ t^'^ otherwise) and T > (^ t^"^ for all i and 
j)r3 A boundary-to-bulk crossover can be observed for sites near the end of a semi- infinite chain at T = oo.a 

The correlation functions {Sf{t)Sj) have a much more complicated structure in the fermion representation. With 
the fermionic identity 

(-1)4'^^ =(4+a,)(4-afc) (2.6) 

applied to ( |2.2| ), this correlation function may be expressed in terms of the auxiliary operators Ak = a| + Ofc and 
Bk = a[, — flfc as follows: 

{Sf{t)S-j) = ^{Ai{t)Bi{t)A2{t)B2{t)...A,^i{t)B,^i{t)Mt)A,BiA2B2...A,^iB,^iA,) . (2.7) 

This expectation value of a product of 2(i + j — 1) fermion operators may be expanded via Wick's theorem in terms 
of more elementary expectation values. The result is most compactly expressed as a Pfaffian:E3 
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\{A,{t)B,{t)) {A^{t)A2{t)) ■■■ (Ai(t)Ai) {A,{t)Bi) 
{Bi{t)A2{t)) ••• (Bi(t)Ai) (Bi(i)Bi) 



{M{t)A,) 



The square of the Pfafhan is equal to the determinant of the antisymmetric matrix with the elements of (2 
the diagonal. The matrix elements can be evaluated from the expressiona23 



{A,it)Ai 



N 



— - sin kj sin fcZ[cos e^t — i sin e^i tanh ^—^] 

k 



(2., 



above 



(2.9) 



{Aj{t)Bi) = — sin kj sin kl [i sin Skt — cos Skt tanh 

k 



(2.10) 



and the relations 



{B,{t)Bi) = -{A,{t)Ai); {B,{t)Ai) = -{A,{t)Bi) 



(2.11) 



All elements ( |2.9D with odd j — I and all elements (2.10) with even j — I vanish. For t = Q the elements ( |2.9| ) are, 
in fact, zero for all j ^ I. All results for the correlation functions {Sf{t)S^) of the XX model presented in the 
following have been derived via numerical evaluation of the determinant associated with (2.8) for systems with up to 
N = 10000 sites. The results are not subject to finite-size effects on the time intervals shown, except where this is 
explicitly stated. 



III. BULK REGIME 



The spin correlation function_/S^(t)S'f^„) of Hxx at T ■ 
identically vanishing for n ^ 0:tj1ij 



oo is surprisingly simple: a pure Gaussian for n = and 



-(5„.oexp 



(T = oo) 



(3.1) 



No physical argument was ever furnished to explain this peculiar decay law. At T = that same spin correlation 
function exhibits a power-law decay, of which the leading term in an asymptotic expansionJ13 



V2A^ 



(r = o) 



(3.2) 



with A = 2^/^^ exp[3C'(— 1)] = 0.64500248..., reflects the Luttinger liquid natumof the ground state of Hxx- Further 
terms of that expansion are known for general n and many more for n — O.EjIlij 

Until recently it was not at all clear whether the Gaussian and power-type decay laws persist at any finite nonzero 
temperatures. Then Its, Izergin, Korepin, and SlavnovEj established on a rigorous basis that {Sf{t)Sf_^.n) decays 
exponentially for < T < oo, i.e. more slowly than ( p.l| ) yet more rapidly than ( |3.2[ ). From the solutions of the 
completely integrable discrete nonlinear Schrodinger model, which is related to Hxx, they were able to derive the 
following expression for the two-spin correlation function: 



exp /(n, 0) , n/Jt>l 
t^"^ ex.pf{n,t) , n/Jt < 1 



(3.3) 



with 
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1 r 

f{n, t) = — / dp\n — Jt smp\ In | tanh 
27r J-^ 



Jcosp 
2T ' 



(3.4) 



1,1 1 

— m tann — — 

27r ' 2T 



(n/Jt)2 



(3.5) 



valid in the space-like {n/Jt > 1) or the time-like {n/Jt < 1) sector of the long-time {Jt — > oo) and/or long-distance 
(n —^ oo) asymptotic regime. The function ,f{n,t) is negative and monotonically decreasing with increasing T; it 
diverges logarithmically at T = oo, thus signalling the change in decay law. 

In the high-temperature regime, the result ( |3.3| ) for the autocorrelation function {n = 0) can be brought into the 
more explicit form 



{smst 



exp 



2Jt. ^ 2T. 
1 + ln — 

TT J 



(1 < T/J < oo). 



(3.6) 



What remains to be filled in for the bulk spin correlation functions (Sf (t)S'f_,_„) is to connect the exponential decay 
to the Gaussian decay (n = 0) or the identically vanishing result (n ^ 0) in the high-temperature limit, and to the 
power-law decay in the low-temperature limit. These connections are realized by crossovers between different decay 
laws at short and long times and can be investigated systematically in the data of finite systems. The salient features 
of the crossovers are described in Figs. and || for autocorrelations (n = 0) and in Fig. ^ for paircorrelations (n 7^ 0). 



A. Autocorrelations 



The modulus-squared spin autocorrelation function \ {Sf{t)Sf)\'^ is plotted logarithmically in Fig. |l| for six values 
of T/J, all in the high-temperature regime (solid lines). The bulk character of these results for site z = 49 of a chain 
with N — 100 spins has been ascertained by comparis on w ith the results of longer chains (with up to iV = 1000 spins). 
The (parabolic) dashed line represents the Gaussian (3.1) - the exact result for T = 00. 

We observe that the Gaussian behavior persists at finite T over some range of short times. That range shrinks with 
decreasing temperature. From the common short-time parabolic shape the individual lines take off like a bundle p£ 
tangents, which represent the exponential character of the long-time decay. The crossover takes place quite suddenly.Ell 
The slight wiggles in the high-temperature data will turn into stronger oscillations in the low-temperature regime as 
we shall see. The observed decay rate in the exponential regime as represented by the slope of the tangent lines 
decreases monotonically as the temperature is lowered. It is well matched in each case by the slope of the adjacent 
dot-dashed line, which represents the decay rate 



1 
r 



2J 

TT 



1 



, 2T' 
In — 
J 



(T/J > 1) 



(3.7) 



of the asymptotic result _([3.6| ). 

The solid lines in Fig. p| show the same quantity as in Fig. |l| but now for three values of T/ J in the low-temperature 
regime. Here the short-time Gaussian behavior has disappeared from the scene. The steepest curve corresponds to 
the highest temperature (2T/J = 1.0). The exponential nature of the decay (with mild oscillations superimposed) 
is now realized even at relatively short times. The average slope is perfectly consistent with the decay rate inferred 
from the asymptotic result (^^), represented by the slope of the adjacent dot-dashed line. 

At lower temperatures a crossover between exponential decay and power-law decay makes its appearance. The 
power-law behavior is first seen at short times. In the center curve of Fig. |, the crossover takes place prior to 
Jt — 20. At longer times the exponential decay is still clearly visible, and the rate of decay agrees well with the 
asymptotic rate 



- = T 

T 



{T/J < 1) 



(3.8) 



extracted from (3^) (see adjacent dot-dashed line). The top curve in Fig. || represent s th e T — result, which has 
been investigated in previous studies. The dashed line shows the asymptotic power law ( ^.2| ) for n — 0, which matches 
the data shown here extremely well (except for the oscillations). 

The onset of finite-size effects at longer times is shown in the inset to Fig. |^ for two cases. The rebound of the 
T = correlation function at Jt ~ 100 can be interpreted as the echo from open ends of the ballistically propagating 
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fermions. In the T > case, the peak at Jt ~ 100 is absent because of destructive interference. The first echo now 
occurs at Jt ~ 200, the time it takes a pulse to move through the system twice. The speed of propa gati on which 



determines the echo time is given by the maximum fermion velocity = max[defe/c?fc]. For the dispersion (2.5) of Hxx 
we have v = J. In the presence ot|anisotropy as realized in the Hxy with ^ Jy, acquires a gap, v decreases, 
and finite-size effects set in later.c3 These echo effects occur at all temperatures. They are easy to recognize in all 
data produced for this study. 



B. Paircorrelations 



In Fig. H we show a logarithmic plot of K^f (i)S'f^„)p for n—A, 9, 14, 19 at the intermediate temperature 2T/J — 1. 
We observe that each function is almost perfectly constant up to a time Jtn — n, where it bends smoothly into 
exponential decay with superimposed oscillations. The decay time does not show any significant dependence on 



n. The inverse decay time predicted by (3.3) for the asymptotic regime of the u ppe rmost curve is given by the 



slope of the dashed line and matches our data very well. A numerical analysis of (3^) shows that the asymptotic 
decay time increases slightly with increasing n at fixed temperature. The linear variation with n of the intercepts at 
Jt = in this logarithmic plot reflects the well-established exponential decay of the equal-time correlation function 

(^f^,%)~exp[-n/e(r)]. 

The inset to Fig. shows again the curve n = 19 of the main plot along with curves for the same correlation function 
at different temperatures. Now the crossover between the space-like and the time-like regime occurs at one common 
value of Jt. In the time-like regime, the slope changes from one curve to the next, which reflects the T-dependence 
of the decay time, while the variable intercept in the space-like regime reflects the T-dependence of the correlation 
length. 

The correlation length r^/T) is known to diverge algeb raically, ^ 1/T, at T = and to vanish logarithmically, 
l/ln(T), at T = ooBa^Ea We have noted in Sec. Ill A that the decay time t(T) also goes to zero logarithmically 



at T = oo and exhibits the same power-law divergence at T = 0. Figure ^ shows both the inverse correlation length 
and the inverse decay time plotted versus temperature. In the XXZ model, the two quantities are expected to have 
more distinct temperature dependences. The numerical results of Rcf. ^ indicate that r(T) stays nonzero at T = oo, 
whereas ^(T) is expected to vanish in that limit as it does in the XX model. 

IV. BOUNDARY EFFECTS 

Here we investigate spin autocorrelation functions {Sf{t)Sf) of a semi-infinite chain for sites i = 1, 2, . . . beginning 



at the boundary. The actual calculations are performed for sites i near one end of a finite open chain (2.3) of N spins. 
However, none of the results presented are affected by the far end of the chain. 

The long-time asymptotic decay of {Sf{t)Sf) in the boundary regime of a semi-infinite chain is different at zero, 
finite nonzero, and infinite temperatures. It is fastest at T = oo and slowest at T = 0, like in the bulk regime, but 
instead of seeing transitions from Gaussian to exponential to power-law decay, we now observe transitions between 
three types of power-law decay. 

At infinite temperature, the power-law long-time asymptotic decay has the form 

(s-nos-f ) - (r = c5o) (4.1) 

with a site-dependent exponent. This result was iniprred from an exact calculation for sites i = 1,2,. ..,5, and 
presumably holds for all sites of a semi-infinite chain.a The poweE-iaw decay at zero temperature as obtained from a 
perturbational treatment of the Pfaffian ( ^.8[ ) is site- independent ,t3 

{Snt)Sf)^t-' (T-O), (4.2) 

and different from the decay in the bulk regime. Our numerical analysis of the same determinantal expression at 
finite nonzero temperatures strongly indicates that the long-time asymptotic decay is a power-law with yet a different 
site-independent exponent: 

{S^it)S^) r^t-"^/^ (G<T<c5o). (4.3) 

An analytic calculation, which confirms this decay law, is presented in Appendix A. The five curves in Fig. |^ 
represent the quantity \{Sf (t)Sf)\'^ for i = 2 at various temperatures in a log-log plot. The data are subject to strong 
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oscillations, which makes it hard to distinguish the different decay laws in this graphical representation. Therefore 
we have smoothed the d ata at Jt > 40 as described in the caption. The three different slopes of the dashed lines 
represent the decay laws (4.1) - (4.3), i.e. for i — 2 (bottom line), t^^f^ (intermediate lines), and (top line). 

The decay laws of {Sf{t)Sf) at a given temperature and for a given site on the semi-infinite chain may undergo 
one or several crossovers. Here we have to deal with bulk-to-boundary crossovers in addition to crossovers between 
different temperature regimes. We first look at the two types individually and then in combinati on. 

At T = oo we ob serv e a bulk-to-boundary crossover, i.e. a crossover from Gaussian decay (3.1) at short times 
to power-law decay (4T) at long times. This is illustrated in Fig. || for the sites i = 2, . . . , 11 near the end of a 
semi-infinite chain. The time Jtc marking the onset of the crossover depends linearly on the distance of the site i 
from the boundary as shown in the inset to Fig. |^. A corresponding bulk-to-boundary crossover between the power 
laws and t-i takes place in {Sf{t)Sf) at T = 0. 

Now we keep the site fixed at « = 2 close to the boundary and vary the temperature. At high-temperatures 
{T/J ^ 1) we can observe an infinite-to-finite-T crossover similar to the one portrayed in Fig. |l| for the bulk regime. 
But here in the boundary regime, it is a crossover between two power laws: and t~^^'^. This crossover is illustrated 

in Fig. 0. The solid lines represent smoothed data at 2T/ J = 10'^ (top) and 2T/J — 10^ (bottom) in a log-log plot. 
It shows how the crossover is shifted to longer times as the temperature approaches infinity. A zero-to-nonzero-T 
crossover between the respective power-laws and t~^/^ can be observed at low temperatures (T/J ^ 1). This is 
illustrated in the inset to Fig. 0. Again the crossover is different from the corresponding (power-law-to-exponential) 
crossover in the bulk regime (Fig. g). 

The two types of crossover in {Sf{t)SJ) of the semi-infinite XX chain, which we have described separately, namely 
the bulk-to-boundary crossover (Fig. P) and the infinite-to-finite-T or zero-to-nonzero-T crossover (Fig. may 
actually occur in one and the same data set. One case in point is demonstrated in Fig. ^ It shows a logarithmic plot 
of \ {Sf{t)Sf)\'^ at high temperatures (T/J ^ 1) for a site of the infinite chain (solid lines) and a site not too far from 
the end of a semi-infinite chain (dashed lines). The solid lines, which represent data already shown in Fig. l|, ex hibi t 
the f ami liar infinite-to-finite-T crossover pertaining to the bulk regime. The same crossover between decay laws (3.1) 
and (3^) is also observa ble i n the dashed lines, but here it is followed by the bulk-to-boundary crossover at finite T, 
i.e. between decay laws (3.6) and ( |4.3| ). 

Note that in the lowest dashed curve, which corresponds to the highest temperature, the intermediate decay law 
has virtually disappeared. H ere the two crossovers overlap. At still higher temperatures their order is reversed. 
The Gaussian decay law ( |3.l| ) crosses over (bulk-to-boundary) to the power- law (4.1), which in turn crosses over 
(infinite-to-finite-T) to the power law (4.3). Corresponding crossover combinations take place in the low-temperature 
regime. 
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APPENDIX A: DECAY LAWS IN THE BOUNDARY REGIME 



The i""^/^ decay law of the spin autocorrelation function {Sf{t)Sf) in the boundary regime at finite nonzero T can 
be reco vere d by t he fol lowing analysis of the Pfaffian (|]^) and its elements. For ^ oo all non-vanishing expectation 
values (2.9) and (2.10) are equal to 



hit) = -i-y-' / dx Vl-x^ (7,-i(x) Ui^iix) e^'^' fiPJx), 



(Al) 



where f{x) = (e^ + 1) ^ is the Fermi function and Un{x) is a Chebyshev polynomial of the second kind. The long-time 
behavior of fji(t) is determined by the singularities of its Fourier transform 



cf>,i{x) = Vl-x^ (7,-1 {x)Ui-i{x) fipjx) 9(1 - x"). 



(A2) 



At nonzero T the only singularities of ( |A2| ) are the square-root cusps \J\ — \x\ 8(1 — |a;|) at x ±1, which yield 
fjiif) ~ e^*'^*t~'^/^. At T = an additional singularity in (j)ji{x) at a; = is generated by the discontinuity in the 
Fermi function, but it contributes only to leading order if both polynomials in (A2) are even, i.e. for odd j and I: 

Instead of the Pfaffi an (P^ ) we study the associated antisymmetric (4j — 2) x (4i — 2) matrix. Its determinant is 
equal to the square of (|2.8|) for i = j. That matrix naturally divides into four (2i — 1) x (2i — 1) blocks: 



M A 
A M' 



The blocks M and M' contain only time-independent elements, fji {0), 
fji{t). Block M turns out to have the following general structure: 



(A3) 



and the block A only time-dependent elements, 



M = 






/l2(0) 
-/l2(0) 




/l4(0) 

-/i4(0) 




V 



/l2(0) 
-/l2(0) 




/23(0) 
-/23(0) 




/25(0) 






/23(0) 
-/23(0) 




/34(0) 
-/34(0) 





/l4(0) 
-/l4(0) 




/34(0) 
-/34(0) 




/45(0) 







-/25(0) 





-/45(0) 





(A4) 



Block M' has a similar distribution of nonzero elements. The leading long-time term in the expansion of (AS) is 
the one which contains the smallest possible number of time-dependent elements, i.e. the largest possible number of 
elements from M an d M '. Given the structure of M, it is not possible to pick more tha n 2i — 2 elements from M 
in the expansion of (A3). The same holds true for M'. Therefore, the lead ing term in (A3) contains exactly two 
time-dependent elements fji{t) from A, which explains the decay laws (4^) and (O) at T = and < T < oo, 
respectively. At T = oo all elements of M and M' vanish, an d the asymptotic long^^^iipe behavior is solely determined 
by block A. This leads to the site-dependent decay law ([4.l[) as explained in Ref. 
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The and decay laws of the spin pair correlation function {Si{t)Sj) in the boundary regime can be explained by 
observing that 4(Sf (t)S|> = {-y~'ih{t)f. 
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FIG. 1. Spin autocorrelation function {Sf{t)Sf) in the bulk regime of Hxx at high temperatures. Plotted is the logarithm 
of |(5f (t)S'f )p evaluated at increments Jdt = 0.4 for site i = 49 in an open chain of 7V=100 spins, for 2T/J = 10, 10^, . . . , 10^ 
(solid lines, top to bottom). The exact result (3.1) for T = oo is shown as dashed line. The slopes of the dot-dashed straight-line 
segments represent the decay rate of (p.fij). 



FIG. 2. Spin autocorrelation function {Sf{t)Sf) in the bulk regime of Hxx at low temperatures. Plotted is the logarithm 
of \{Si {t)Si)\ evaluated at increments Jdt = 0.4 for site i = 49 in an open chain of = 100 spins at 2T/J — 1.0, 0.1 and for 
site i = 249 in a system of size A'' = 500 at T = (solid lines, bottom to top). The dashed line represents (3.2). The slopes of 
the dot-dashed straight-line segments represent the decay rate inferred from the asymptotic expression (3.3). The inset shows 
the results at 2T/J = (upper curve) and 0.1 (lower curve) again over a longer time interval and (now in both cases) for 
N = 100 in order to illustrate the onset of finite-size effects in the form of echos due to ballistically propagating pulses reflected 
at the open ends of the chain. 



FIG. 3. Spin paircorrelation function {Sf{t)Sj) in the bulk regime of Hxx at 2T/J = 1, for j = 49 and i — 30,35,40, 
and 45 (solid lines, top to bottom) in an open chain with A'^ — 100 spins. Plotted is the logarithm of \{Si {t)Sj)\'^ evaluated 
at increments Jdt = 0.4 The inset shows the same function in the same representation, but here the sites are kept fixed 
{i = 30, j = 49) and the temperature is varied {2T/J — 0.1, 1 , 10) . The slopes of the dot-dashed straight-line segments 
represent the decay rate inferred from the asymptotic expression (|3.3|). 



FIG. 4. Temperature dependence of the inverse decay time 1/r as derived from the asymptotic expression (3.3) (solid line) 
and from our numerical analysis of large but finite chains (circles). The numerical data were derived from an exponential fit 
to \{Sf{t)Sf)f for i = 49 in an open chain of = 100 spins, evaluated for 20 < Jt < 60 at increments Jdt = 0.4. The dashed 
lines represent the approximations (3.7) and (3.5) to (3.3) for high and low T, respectively. For comparison we also show 
(dot-dashed) the temperature dependence of the inverse correlation length 1/^ as given by (3.J). Data read ofi' Fig. 2 of Ref. 



coincide with the curve shown, as do data derived from our own numerical results for (Sf Sj) / (Sf^iSj) with 
20 in an open chain of A'^ = 100 spins. 
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FIG. 5. Spin autocorrelation function {Sf{t)Sf) in the boundary regime of Hxx at 2T/J — 0, 0.1, 1, 10, oo (solid lines, top 
to bottom). Shown is the quantity \{Sf {t)Sf)\'^ evaluated at increments Jdt = 0.4 (Jt < 40), Jdt = 1 (40 < Jt < 100), and 
Jdt — 10 (100 < Jt < 2000) in a doubly logarithmic plot. The data pertain to the site j = 2 in an open chain of A'^ = 10000 
spins. At times Jt > 40 most of the data were smoothed by the following two-step procedure: (i) retain all local maxima of 
the data set; (ii) eliminate all local minima from the remaining data set. The 2T/J — 0.1 data did not require smoothing. The 
dashed lines represent the power laws t~^ (top), t~^^^ (center), and t~^^^ (bottom). 



FIG. 6. Logarithmic plot of the square of the spin autocorrelation function (Sf (t)S f) in the boundary regime of Hxx at 
T — oo for j = 2, . . . , 11 (top to bottom). The dashed line represents the Gaussian ( p.l[ ) pertaining to the bulk limit {i — » oo). 
The data for i < 5 are obtained from the exact expressions given in Ref. ^ and the data for i > 5 from an open chain with 
N = 100 spins. The inset shows the time Jtc at which the relative deviation of (Sf (t)5'f )^ from the Gaussian first exceeds one 
percent. 



FIG. 7. The main plot shows smoothed data of the same quantity as in Fig. H at two very high temperatures 
(2T/J — 10"^, 10^) and in the inset for two very low temperatures 2T/J = 2 • 10~^, 10~ ). The dashed lines in the main 
plot represent the power laws t~^^^ (top) and t~^^^ (bottom). The dashed lines in the inset represent data for temperatures 
2T/J — 0,0.1 (i.e. the two uppermost solid lines of Fig. 5). 



FIG. 8. Spin autocorrelation function {Sf{t)Sf) in the bulk regime and in the boundary regime of Hxx at high temperatures. 
Plotted is the logarithm of \{Sf {t)Si)\'^ evaluated at increments Jdt = 0.4 for site i = 49 (solid lines) and site j = 11 (dashed 
lines ) in an open chain of TV = 100 spins for 2T/J = 1, 10, 10'^, and 10^. The dot-dashed line represents the Gaussian decay 
law jill) . 
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